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£2 ! Abstract 

O 



In this paper various predictions for the scaling exponents of the Nonlocal Kardar-Parisi-Zhang 



(NKPZ) equation are discussed. I use the Self-Consistent Expansion (SCE), and obtain results 
that are quite different from result obtained in the past, using Dynamic Renormalization Group 

o, 

analysis (DRG), a Scaling Approach (SA) and a self-consistent Mode Coupling approach (MC). 



It is shown that the results obtained using SCE recover an exact result for a subfamily of the 
NKPZ models in one dimension, while all the other methods fail to do so. It is also shown that 
the SCE result is the only one that is compatible with simple observations on the dependence of 
the dynamic exponent z in the NKPZ model on the exponent p characterizing the decay of the 
nonlinear interaction. The reasons for the failure of other methods to deal with NKPZ are also 



discussed. 
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I. INTRODUCTION 

The field of disorderly surface growth has received much attention during the last two 
decades. A multitude of different phenomena such as fluid flow in porous media, propagation 
of flame fronts, flux lines in superconductors not to mention deposition processes, bacterial 
growth and "DNA walk" |l[ are all said to be related to the famous Kardar-Parisi- Zhang 
(KPZ) equation 2]. Therefore, further understanding of the behavior of the KPZ equation 
has a broad interest in the fields of nonequilibrium dynamics and in disordered systems. 

The KPZ equation for a growing surface is 

dh(f.t) „ 9 , A.„, n2 ,_ N 

K dt J = vV 2 h + - (V/i) 2 + V (f, t) , (1) 

This equation describes the height fluctuations of a d-dimensional interface, where h (r, t) is 
the height at r measured relative to its spatial average, v is the diffusion constant and r\ (r, t) 
is the fluctuation of the rate of deposition. The noise 77 (r, t) has a zero mean ((77 (r, £)) = 0) 
and it satisfies 

(7/ (f, t) 77 (7, 0) = 2D S d (r- r 1 ) 5 (t - t') , (2) 

(white noise) where d is the substrate dimension and D specifies the noise amplitude. 

Together with the great success of the KPZ equation to describe many growth models and 
phenomena, there has been a growing pool of data that is not well described by KPZ, and 
triggered further research. One of the first classes that belongs to this non KPZ behavior is 
the well known Molecular-Beam- Epitaxy (MBE) class (sometimes also called the conserved 
KPZ equation) J3|-J9(. This class is distinct from KPZ in having surface diffusion as the basic 
relaxation mechanism. However, the modified behavior introduced by the MBE equation is 
not at all sufficient to account for all the rich non KPZ experimental data in the field. 

A different line of research suggested that the basic growth equation should not be 
changed. Instead, the white noise that appears in the original KPZ equation should be 
correlated - either temporally or spatially [10]-[14]. This approach was indeed quite success- 
ful. However, it still failed to give a good account for all the measured scaling exponents. 

Recently, some researchers suggested that incorporating the long-range nature of inter- 
actions in the gro wing surface is necessary for a proper description of many systems such 
as colloid systems [15j . paper-burning experiments |li| or protein deposition kinetics [17] . 
Following this basic intuition Mukherji and Bhattacharjee [1S| developed a Langevin-type 



equation with a nonlocal non-linearity, thus going beyond the KPZ local description of inter- 
actions. They studied originally the white noise case that was later generalized to spatially 
correlated noise by A. Kr. Chattapohadhyay |L9j. To be more specific, the equation they 
studied was 

Oh (t t) I r 

— ±2-L = u\/ 2 h(f,t) + - dr'g(r')Vh(f+r',t) ■Vh(f-r',t) +r](r,t), (3) 

where the kernel g (r) represents the long-range interactions. They take g (r) with a short- 
range part \o5 d (r) and a long-range part ~ X p r p ~ d , or more precisely in Fourier space, 
g (q) = A + A p g~ p . The noise has again a zero mean, but it is allowed to have spatially long 
range correlations, characterized by its second moment 

(rj (r, t) v if, t)) = 2D, \r- ff^ 5 (t - t') , (4) 

where the case of white noise is restored in the limit a = |lj| . 

In the next section I will present current theoretical predictions for the scaling exponents 
of the NKPZ model. It will become evident that all the previous results are inconsistent with 
an exact one-dimensional result obtained in the past. Then, in section III the Self- Consistent 
Expansion (SCE) approach will be applied to this model. Interestingly, this method yields 
results that are consistent with the exact one- dimensional result. Then, results for higher 
dimensions as well as for general p's are presented. Eventually, a discussion of the reasons 
for the failure of the various theoretical methods is given in section IV, and conclusions are 
drawn in section V. 

II. PREVIOUS RESULTS 

Both papers |l8|-[19( investigated this problem using Dynamic Renormalization Group 

(DRG) analysis, and derived a very complex zoo of universality classes (ref. [19] also applied 

a mode-coupling approach confirming some of the results of DRG). Let us not get into 

all those details, but rather focus on the strong coupling solution (in the KPZ sense |l|) 

suggested by both papers, namely 

_ ( d -2-2p)(d-2-3p) 
ZDRG ~ 2 + (3 + 2-P)d-6-9p ' (5) 

where z is the dynamic exponent. To complete the picture, the roughness exponent, a, can 

be calculated using the modified scaling relation a + z = 2 — p (which actually comes from 



the famous Galilean invariance). Notice that this result does not depend explicitly on a as 
long as a > 0. 

On the basis of the DRG result (eq. EJ) Mukherji and Bhattacharjee tried to explain some 
non KPZ experimental results in one dimension, namely a roughness exponent of a = 0.71 
given in refs. (la|-(la| that implies, according to their DRG calculation, p = —0.12, and 
requires A = 0. As already appreciated by the authors, this result raises some doubts on 
physical grounds, whether a true Aq = is seen in experiment, and where should a negative 
p, which implies anticorrelations 20] rather than correlations that cause the interaction, 
come from. 

This finding is obviously disturbing because if DRG yields accurate results then the NKPZ 
equation cannot be relevant to the physical processes discussed above. As will be shown 
later, the DRG result is far from being accurate, so that the above argument against the 
relevance of the NKPZ equation is not valid. The relevance of the NKPZ equation has been 
doubted, however, for other reasons too, but this will not be discussed here. 

The equation is interesting by itself and generalizes the traditional KPZ prototype of 
nonlinear stochastic field equations that are so abundant in the description of natural phe- 
nomena. Therefore, reliable methods of solution are of great importance. DRG is such a 
general method, but unfortunately the DRG results look suspicious for three reasons: 

First, the expression for the dynamic exponent (eq. (jSJ)) reduces to the well-known result 
of Medina et. al. [IQj z = 2 + ^_ 6 in the limit of p = (i.e. the limit of local KPZ) - a 
result that by itself is meaningful only for the case of d = 1. It is a well-known fact that 
DRG cannot give the strong coupling solution of the local KPZ equation for d > 1 22]. 
This may suggest that the DRG result is not correct also for the nonlocal case, especially 



for d > 1. Furthermore, the mode coupling approach that is 
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mown to be more appropriate 
23j | yielded a strong coupling 



than DRG when dealing with the standard KPZ problem 
result that is different from eq. (jHJ) |l9j| . Yet, it should be mentioned that the mode coupling 
approach |l^| (at least using the simplest ansatz) was not able to produce a non power 
counting solution, as obtained in the local KPZ case [23j ]. 

Second, it is easily seen that the explicit expression for the dynamical exponent (given in 
eq. (0) above) has a singularity for p's that satisfy (3 + 2 _Ps ) d — 6 — 9p s = (for example 
in one dimensionp s = —0.205). It is not clear why should such a singularity appear at all 
and if it appears, then why should it take this specific value. 



Third, from general considerations one should expect the following scenario: as p becomes 
more and more negative, the interaction that couples the gradients is enormously reduced 
with distance (remember that the long-range part of the kernel scales as g (r) ~ \ p r p ~ d ). 
Therefore, one should expect that in this region, the nonlinear term becomes irrelevant, 
and the scaling exponents of the linear theory (also known as the Edwards- Wilkinson (EW) 



theory 



21| ) should be recovered. On the other side, as p becomes larger and larger, the 



interaction becomes more and more relevant, and one might expect that as a consequence 

faster relaxations should appear in the system (faster relaxations imply smaller values of 

the dynamical exponent z). Therefore, it is reasonable that a non-increasing behavior of the 

dynamical exponent (z) as a function of p will be seen. It is not difficult to see that the 

expression for the dynamical exponent (given by eq. ©) does not follow this reasoning. 

I 
Lately, an exact result for the NKPZ model was found [2J|. It turns out that the Fokker- 

Planck equation associated with the Langevin-form of the NKPZ model (i.e eq. above) 

can be solved exactly for a specific sub-family of models in one dimension. More specifically, 

a Gaussian steady-state solution was found for the case Ao = 0, p = 2a when d — 1. This 

exact solution yields the following dynamic exponent 

_ 3-3p 

Zexact c. ■ \y) 

It is easy to see that this result reduces to the well know local-KPZ (that corresponds to 
p = ) result, z = 3/2, in one dimension. 

This exact result is not compatible with the DRG result - either quantitatively nor quali- 
tatively. As mentioned above, DRG is usually considered relatively reliable in one dimension. 
However, this exact result indicates DRG's shortcoming already in one dimension. The in- 
evitable conclusion is that in order to gain insight into the behavior of the NKPZ model in 
any dimension, and for any value of p (the nonlocal parameter) and a (the long-range noise 
parameter) it is necessary to employ more reliable methods. 

Actually, two methods that proved useful in the context of the local-KPZ problem were 
applied recently to the NKPZ problem as well. The first method, is a Flory-type scaling 
approach that was originally proposed by Hentschel and Family [25J] and generalized lately 

I I 

by Tang and Ma [26[ to the nonlocal case. The strong-coupling dynamical exponent obtained 

by using this Scaling Approach (SA) is 

(2-p)(2 + d-2a) 
ZSA = 3 + d-2a • (?) 
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The second method applied to the NKPZ equation is an improved version of the mode- 
coupling approach. More precisely, it is a solution of the mode-coupling equations using 
an improved ansatz that describes correctly the asymptotic behavior of the solution. This 
approach was first proposed by Colaiori and Moore |27l |. and applied by Hu and Tang [2j| 
to the nonlocal case. The strong-coupling exponents using this method are obtained by a 
numerical solution of the following set of transcendental equations 

PAS' d _ d(2-z) 2 2 p 
B ~ nz 2 BI(B,z,p) 

PAS' d p(p-z) 2( 2 ^)/^ 

B 2 " nZ 2 V(^)/P ( } 

where S' d = f£ sin d ~ 2 9d0, T (u) is Euler's gamma function, P = d + 4 — 2z — 2p, B — 2(2 1 _ z - ) 
and I (B, z, p) = f£° (1 — p — 2s 2 ) s 2z_3 exp (—B^' z s^ — s 2 ) ds. Such a numerical solution 
for z as a function of p for d — 1, 2, 3 is presented in FIG [D of ref. [28j |. 

While these two results look reasonable, in the sense that both do not develop any 
singularity as a function of p, and at least the first one (given by eq. (J7J)) describes a non- 
increasing dynamical exponent, it turns out that both results are still incompatible with the 
exact result of z = (3 — 3p)/2 (given by eq. ©)• In addition, the mode-coupling solution 
has a peak for d > 1 meaning that the dynamic exponent z is not monotonous as a function 
of p that is expected from general considerations. Therefore, there is still a need for a reliable 
method to tackle this problem. The reasons given above motivated the analysis presented 
in this paper. 

III. THE SELF CONSISTENT EXPANSION 

In this paper I apply a method developed by Schwartz and Edwards [lj, |29|, |30| (also 
known as the Self-Consistent-Expansion (SCE) approach). This method has been previously 
applied successfully to the KPZ equation. The method gained much credit by being able 
to give a sensible prediction for the KPZ critical exponents in the strong coupling phase 
for d > 1, where, as previously mentioned, many Renormalization-Group (RG) approaches 
failed (as well as DRG of course). It also produces the exact one-dimensional result. I will 
show that in contrast to DRG, the experimental result a = 0.71 is accounted for, when using 
the SCE, by a positive p. Furthermore, I will show that the SCE method yields the exact 



result, for the specific sub- family of models that can be solved exactly in one dimension - 
where other methods fail. In addition, solutions for any d, p and a are obtained. These 
solutions are consistent with the expected qualitative scenario presented above. 

As will be seen immediately, another remarkable advantage of the SCE method is the 
minor changes needed in order to generalize the result of local KPZ with uncorrelated noise to 
include nonlocal interactions as well as spatially correlations in the noise. The above implies 
a second important motivation for this paper, namely a demonstration of the robustness of 
the SCE method as well as its mathematical coherence and consistency. 

The SCE method is based on going over from the Fourier transform of the equation in 
Langevin form to a Fokker-Planck form and constructing a self-consistent expansion of the 
distribution of the field concerned. 

The expansion is formulated in terms of 4> q and u> q , where <f> q is the two-point function in 
momentum space, defined by <p q = (h q h_ q ) s , (the subscript S denotes steady state averag- 
ing), and u q is the characteristic frequency associated with each mode. 

It is generally expected that for small enough q, <p q and uo q are power laws in q, 

<P q = Aq~ T (9) 

and 

u q = Bq\ (10) 

where z is just the dynamic exponent, and the exponent V is related to the roughness 
exponent a by 

V — d , „ 

a = -2~. (11) 

The method produces, to second order in this expansion, two nonlinear coupled integral 
equations in <p q and u q , that can be solved exactly in the asymptotic limit to yield the 
required scaling exponents governing the steady state behavior and the time evolution. 

In the following I will find the critical exponents of the NKPZ equation with spatially 
correlated noise. I consider the NKPZ equation (eq. (|HJ)) in Fourier components 

dh 1 

-^ = -u q hq j= Y^ M qem h e h m + r) q , (12) 



where Q is the volume of the system, to be taken eventually to infinity and r\ q is the noise 
term that satisfies: 

M*)> = 0, (13) 

(Vq (t) V q > (*')) = 2D q S d (q + <?) 5 (t - t') , (14) 

and M q i m , v q and D q are defined by 

A, 



M„ 



~qlm = —]=<l P (/■ m) Sq,e+r. 



v q = vq 



D q = D q 



-2a 



(15) 



Thus, the NKPZ equation is of the general form discussed in refs. [29, 30], where the SCE 
is derived. Working to second order in the expansion, this yields the two coupled non-linear 
equations 



7-. j . n v^ M q e m M q £ m (j)£(fi m ^-^ M q i m Mi mq (j) m (j) q 
D q -v q( f> q + 2^ ,,,,,, 2 L 



Ira U 1 ~^~ Ue + U ' n 



i.m 



UJ q + Ui + U r , 



o v-^ M q n m M mtq (j)i(f) q _ 

Ira ^q + U £ + U m 



and 



Vq ~Uq 



2E M « 

£,m 



Mt mq (j) m + M mtq 4>t 



0, 



(17) 



where in deriving the last equation I have used the Herring consistency equation 31]. In fact 
Herring's definition of ., is one of many possibifities, eacfi leading to a differ** cogency 

equation. But it can be shown, as previously done in 30], that this does not affect the 
exponents (universality). 

A detailed solution of equations (fTH|) and (|T7jl in the limit of small q's (i.e. large scales) 
is performed (as in refs. [14, 30]), and yields a very rich family of solutions for any d, p 
and a. These solutions are given in Table |U However, before getting into all the details 
let me focus on the strong coupling non-power counting solution obtained by SCE. The 
reason for focusing on this solution is that it is interesting to compare it with the solutions 
obtained by other methods, and with the experimental results. This solution is determined 
from the combination of the scaling relation z = + ~ ~ p , obtained from eq. (j!7j) . and the 
transcendental equationF (r, z, p) = 0, where F is given by 



F(T,z,p)= - d d t 



t- (e - t) 

t z + e-i\ Z + 1 



(«•*) 



e-t 



t\-P 



t~ l +e 



(M 



+ 



+ / d d t- 



e-t 
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t z + 



t 



-r 



+ 1 



t 



(18) 



and e is a unit vector in an arbitrary direction. 

This solution is valid as long as it satisfies the following conditions: d < 3r (d, p) — 4 + 
min {2p — 4a, — 2p} and d < T (d, p) — 2 |p|, where r (d, p) is the numerical solution of the 
equation F (T, z, p) = - when such a solution exists (how to obtain the conditions above 
see ref. [30]). 

It turns out that for d — 1 the equation F (T, z (r) , p) = is exactly solvable for any p, 
and yields T = 2 + p and z = -^rp (it can be checked immediately by direct substitution). 
In this case the three validity conditions read — 1/3 < p < 1 and p > ^p-. By using eq. 
(jlljl I translate the results into a = ^- and z = ^-^-- It is straightforward to check that 
this solution reduces to the exact KPZ results of a = 1/2 and z = 3/2 in the limit of p = 0. 
Moreover, this solution recovers the exact solution obtained in ref. [24], and by doing that 
clearly does better than all the other results presented above (eqs. ©,© and (jEJ)). 

It should also be mentioned that for d > 2 such an exact solution for F (T, z (T) , p) = 
as a closed analytical expression cannot be found, and one has to solve numerically the 
equation in order to find T (d, p) (such results are presented below). 

As can be appreciated, the full description of the results given in TABLE 1 is quite rich 
and may look confusing. In order to gain more insight into this seemingly long list of possible 
phases I concentrate on the case of white noise (a = 0), and I describe all the possible phases 
(as a function of p) for d — 1,2, 3. 

In one dimension, the following scaling exponents are obtained: 



5-2p 

3 
3-3p 

2 



2p 



p<-\ 

-i<p<- 

-|<p<l 

p>\ 



and 



p<-\ 



5-2p 
3 

2 + p 

2 



"2 < P < — 5 
-\<P<1 
P >\ 



(19) 



It can be seen that the dynamic exponent obtained using SCE is continuous, has no 
singularities and is non-increasing as a function of p. In addition, the intuition that the 
scaling exponents of the linear theory should appear in the limit of very small p's is verified 
and specified, namely for p < — | the Edwards- Wilkinson (EW) exponents of z = 2 and 
a = 1/2 are obtained. 



An important observation regarding the one- dimensional results is that for every region 
of p's, well defined values of z and T are obtained, so that only one phase is possible when 
p is specified, and no phase transition as a function of the dimensionless coupling constant 
(that is defined by g 2 = ^-^) exists. 

In two dimensions, the following scaling exponents are obtained: 



2 p< 

2 - 2p p > and T 

^f^ + 2-p A< P <pl 



2 p< 

2 p>0 , (20) 

r (2, P ) p^ 2 <p<^ 



where P2 an d P2 are the lower and upper bounds on p (respectively) for which a solution 
for the transcendental equation F (T,z (T) , p) = yields a dynamic exponent that does not 
exceed z = 2 and z = 2 — 2p. To be more specific, p f 2 = —0.395 and p% = 0.305 (the dynamic 
exponent is also presented in FIG. 1 (b)). 

An important observation regarding the two-dimensional result is that for values of p 
within the region pi < P < P2 there are two possible solutions (except for p = 0). This can 
be easily seen in FIG. 1 (b) - where two branches appear in that region. This means that 
a phase transition is possible, as a function of the strength of the dimensionless coupling 
constant in the theory, between a weak coupling solution (the upper branch in the figure) 
and a strong coupling solution (the lower branch in the figure). 

This phenomenon is surprising in view of the known fact that d = 2 is the lower critical 
dimension of the local KPZ theory. This means that for dimensions higher than 2 a phase 
transition becomes possible in the local problem, but in 2D is only the strong coupling 
phase. This fact is indeed respected by the results presented above, because for p = 
(that corresponds to the local theory) such a phase transition is not possible. However, 
surprisingly if one takes either positive or negative p's (even slightly higher or lower than 
zero) such a phase transition becomes possible. This finding is explained by the fact that 
the lower critical dimension is lowered in the presence of long-range interactions. The key 
point is that for positive p's (p > 0) the relevant "weak coupling" phase is a phase with a 
dynamic exponent of z = 2 — 2p, and this phase becomes possible for d > 2 — p. In addition, 
for negative p's (p < 0), a different "weak coupling" solution is possible - this one with a 
dynamic exponent of z = 2, and this phase becomes possible for d > 2 + 2p = 2 — 2|p|. 
Therefore, once we encounter a non-zero p the lower critical dimension is reduced (for two 
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different reasons), and one of the weak coupling solutions becomes possible. In that case, a 
phase transition can also take place. 

In more than two dimensions, the scaling exponents follow the scheme: 



2 p< 

2 - 2p p > and T 

^^ + 2-p Pd <p<p u d 



2 p< 

2 p>0 , (21) 

, ^ y y d ^ r ^ y d V T (d, p) p e d < p< p u d 

whenever a solution for F (T,z (T) , p) = exists. Notice a difference from the two- 
dimensional case in that the weak-coupling exponents are also possible when p = 0. I 
have performed the detailed calculation for d = 3, and it is shown in FIG. 1 (c) (where 
Pg = —0.016 and pg = 0.246). It can be seen that no big quantitative difference from the 
two-dimensional case is present. 

IV. DISCUSSION 

The results I obtained using the self-consistent expansion (SCE) are obviously different 
from those obtained by all the other methods. First, the SCE is able to reproduce the 
exact one-dimensional result, where all the other methods fail. Second, the SCE results are 
compatible with the expected behavior of the dynamic exponent for p — > — oo (i.e. recovery 
of the EW exponents) as well as with the expected non- increasing behavior of the dynamic 
exponent as a function of p. None of the other results mentioned before recover the EW 
behavior for large negative p's, and only one such result shows a non-increasin g b ehavior of 
the dynamic exponent as a function of p - namely the result presented in ref. J26J]. 

At this point, it might be interesting to return to the experimental result of refs. |l5j]-|l6| 
that were discussed previously, and to discuss them in view of the new results derived here. 
It turns out that the SCE method yields p = 0.42 (p = 2a — 1) for the case of a = 0.71 
(at least if the noise is white, or more precisely if the exponent that describes the decay of 
spatial correlations a is not large). This result is physically more reasonable than the result 
p = —0.12 suggested in ref. [18j . In addition, since p > 0, it does not require to impose the 
strict requirement A = 0, that is also problematic physically. 

It remains a mystery why all the methods that were mentioned at the beginning of the 
paper failed to reproduce the exact one-dimensional result, while they all yield the exact 
result for the local KPZ case (p = 0). The answer to this puzzle lies, in my opinion, in the 
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fact that all those methods fail to notice that the nonlocal nonlinearity generates a nonlocal 
relaxation term under renormalization. This means that once a nonlocal nonlinearity is 
introduced into the equation, effectively a certain fractional Laplacian is present as well. 
This fact can be observed when applying the SCE method, in the "one-step renormalization 
stage" that is part of the solution (see for example eqs. (6.8)-(6.10) in ref. |2I|). Then, 
in addition to the g 2 -relaxation term a g 2_2p -relaxation term appears, and it obviously 
dominates the dynamics in the large-scale limit (i.e. q — > limit) when p > 0. I think 
that the failure to realize this fact prevented the other methods from recovering the exact 
one-dimensional result. 

The last statement can be checked for the methods mentioned above. My prediction 
is that if one adds by hand the correct relaxational term, namely q 2 ~ 2p , then at least in 
one dimension the right answer will be obtained for the nonlocal model as well. This idea 
suggests another important advantage of SCE, since a straightforward application of SCE 
yields the extra fractional Laplacian, and it does not have to be thought of and added 
beforehand. 

V. CONCLUSIONS 

In this paper several results for the nonlocal Kardar-Parisi-Zhang equation have been 
discussed. It has been shown that those theoretical predictions for the scaling exponents of 
the NKPZ equation are inconsistent with a recent one-dimensional exact result [2J| that is 
possible for that model. I also claimed that these predictions are not compatible with the 
expected behavior for the dynamic exponent z as a function of the nonlocal parameter p 
from general considerations presented above. 

Then, in order to derive reliable results that are consistent with the exact one-dimensional 
result an alternative method, namely the self-consistent expansion (SCE), has been used. 
The results obtained using SCE were also consistent with the expected p-dependence of 
the dynamic exponent. In addition, the SCE is able to give sensible predictions for higher 
dimensions, as well as for any value of p and a. All the possible phases found using SCE 
are summarized in Table I. 

I then discuss the implication of this calculation to certain experimental results |15J]-|1(t| 
that might be well described by a NKPZ model. At the end I also discuss possible reasons 
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for the failure of all the other methods to recover the exact one-dimensional result, and 
to yield sensible predictions, namely the generation of a fraction Laplacian term in the 
equation under renormalization. This observation may trigger future work on the "fixed" 
NKPZ model that includes such a term from the beginning, using methods such as mode- 
coupling. It would then be very interesting to compare between the various results in higher 
dimensions. 

This situation, along with previous results obtained in 1J, |29|, |3(|, suggest that the 



SCE method is generally a successful method when dealing with such non-linear Langevin 
equations. 

Acknowledgement: I would like to thank Moshe Schwartz for useful discussions. 
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z 


r 


Validity 


2-2/9 


2-2p + 2a 


a > p> and d > 2 + 2a - 3p 


2-2/9 


2 


p > , a < p and d > 2 — p) 


2 


2 + 2p 


a < p < and <i > 2 + 4/9 


2 


2 + 2a 


p < , a > p and one of the following: 

1. d > max{2 + 2o- + /9, 2 + 4cr}. 

2. o- > \p and d > 2 + 2a + 2p. 

3. a <\p and 2 + 4ct < d < 2 + 2a + p. 

4. 2 + 2a + 2p < d < 2 + 4a. 


d+4-2p-2cr 
3 


d+4-2p+4a 
3 


a > p , d < 2 + 2a - p - 3\p\ and d + 4 - 2p + 4a > 3F (d, p) 


d+4-r (d,p)-2p 

2 


Sol. oiF(F,z,p) = 0, 
denoted by Tq (d, p) 


d < 3r (d, p) - 4 + min {2/9 - 4cr, -2p} 
and (i < Tq (d, p) — 2 \p\ 



TABLE I: A complete description of all the possible phases of the NKPZ problem, for any value 
of d, p and a. The first two columns give the scaling exponents z and T for a particular phase, 
and the third column states each phase's validity condition. Note that Tq (d, p) is the numerical 
solution of the transcendental equation F (T, z, p) = with the scaling relation z = 2 ~ P ■ 
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The nonlocal parameter p 




The nonlocal parameter p 




The nonlocal parameter p 



FIG. 1: The values of the dynamic exponent z as a function of the nonlocal parameter p, for 
uncorrelated noise (a = 0). (a), (b) and (c) are for d = 1, 2 and 3 respectively. 
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